setwd("D:/杨熊冰/R语言/实验报告/data for class3")
#5.1 数据分布
ADdata <- read.table("ADdata.txt") #读取实践3存储的文件
GSEdata <- read.table("GSE67835_data.txt") #读取实践3存储的文件
ADdata <- as.matrix(ADdata) #转化成矩阵
hist(ADdata) #直方图
hist(ADdata[ADdata != 0]) #去除0之后的数据
hist(log2(ADdata+1)) #对数化 +1 是让0可以取对数
hist(log2(ADdata[ADdata != 0])) #去除0之后的数据对数化

GSEdata <- as.matrix(GSEdata) #转化成矩阵
hist(GSEdata) #直方图
hist(GSEdata[GSEdata != 0]) #去除0之后的数据
hist(log2(GSEdata+1)) #对数化 +1 是让0可以取对数
hist(log2(GSEdata[GSEdata != 0])) #去除0之后的数据对数化
hist(log2(GSEdata[GSEdata > 5])) #>5的数据对数化
